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Contact with 
has indicated the 
pressing, yet the 
concerning the pr 



cone 1'achine Simulation llodels (war games) 
t the need for "random numbers" is real and 
information available to the war gamers 
operties of these numbers is often quite 



vague and limited. 

The average of the results of several runs of a war game 
is often used to evaluate the many relationships tested by 
the game. Since random numbers are the device that determine 
the yes/no answer at the many decision poirits of the game, 
the faithfulness with which these numbers exhibit the proper- 
ties desired by the war gamer directly affect the validity of 
the results. 



In view of the critical role in war gaming played by 
these numbers it seems useful to provide a tested method of 
producing the numbers, and to provide also the means of test- 
ing any other number-producing methods considered useful. 

Random numbers have an important application in statis- 



tical sampling, war gaming, and other monte carlo techniques. 
The applications mentioned in this paper will be primarily in 



:he war game area. 

This thesis was written at the United States have. 1 Tost- 



durins the Deriod 



graduate School, Uonterey, California, 



January-' Vo, , 1'*'?. I tx i 

for for his encourage 
noting a s faculty a ^ 7 * ' or, 



ndeht or’ to * rofsssor .’ac 1 ' r . . 

. er.t and capable guidance while 
a 't 1 to Professor Thoma s h 



Cho^vc’ for hi s ve 1 vable acrirtance <.-s second reader. I 
wird’ to tharh t’w> '-’aval lost graduate dchool Conoutev 'fencer 
staff for their patience and efficiency in accepting and 
process 5 tv* a multitude of programs compiled during the writ 



in? of this thesis. 

This thesis was ’written irimari by for the reader with 
had: ground in calc ll is and statistics, but may be useful to 
anyone with a basic ’-novledge of the use of computers for 



war 



The U. T . -;:vy has procure ' 1 high s^eed digital computers 
and ad .. pted automatic data processing AVN) to nc-ct the t'e- 
naru's of the increasing complexities • e modern warfare and 
weapon r«> stems . In particular, such computers and i’ll find 
ay plication in the Naval Tactical Tate ‘"hasten and in Command 
an - 1 Control Tenters which h, ve been set up to provide -■ effec- 
tive direction of military forces. These military activities, 
with automatic data processing capability, also have a poten- 
tial capability for war pa min; an 1 machine r. i uvula ti on with 
the corvputers . ! Ian}’ activities are engaged in the writing, 

playing an .' 1 evaluation of war awes at this time. It seems 
lively that war .gaming will soon become or.e of the routine 
functions of the commands with A T V rs capability. 

To assist in meeting the educational requirements of this 
situation, formal instruction in war gaming u-ing digital com- 
puters is under development at the U.C. Naval Academy and t l, e 
U. ^ . Naval Postgraduate "chool. Fleet Fropramrning Center:, 
have been established for direct support of activities employ- 
ing automatic data processing and conducting machine simula- 
tion war games. 

Through the formal instruction in war gaming at the N.~. 



Nas al “hostgre hate School it he; a become ac or. rent t 1-r t 









are none initial problem" in ^hc study of war gaming. One 
problem lies ir. I aving precise terminology and notation to 
explain concepts that are largely probabilistic in nature. 
Another lies in the stochastic structure of the game, i.e., 
the decision points in the game which require drawing a ran- 
dom number to determine the occurrence or non-occurrence of 
an event. This provides the requirement for a random number 
generator in the game. A third problem is, given a random 
number generator, are the characteristics of the numbers pro- 
duced by this generator such that the results of the game are 
valid? 

The second, and third problem areas mentioned, that is, 
generation of random numbers, and testing of random numbers 
are the subject of the investigation in this thesis. The 
investigation included a search of the literature and accu- 
mulation of reference material, writing CDC 1604 FCTTTAN pro- 
grams for generation of random normal, exponential and chi- 
square numbers, writing programs for tests of random numbers, 
and running the programs for the generators and tests on the 
Tt*C 1604 computer. 

Sections 1 and 2 of the thesis contain discussions of 
generation of random numbers and testing of random numbers, 
respectively. Section 3 contains the results of tests per- 
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formed on two normal random number generators, an exponen- 
tial and a chi-square random number generator. Section 4 
presents the conclusions reached In this paper. Appendices 
A, B and C provide justification for the methods of genera- 
tion of random numbers and Appendix D reproduces the pro- 
grams in CDC 1604 FORTRAN language for the four generators. 

It is suggested that the reader familiarize himself with 
the explanation of notation and definitions of terms immediate- 
ly following this introduction. It is also suggested that the 
U.S. Naval Postgraduate School thesis "R.andom Number Genera- 



tion on Digital Computers by J. Barron 



be consulted 



for information concerning the generation of uniform random 



numbers . 



DEFINITIONS OF TERMS AND NOTATION 
In this paper a random variable, X, and its associated 
cumulative distribution function F, will be identified as 
(X,F) or, in some instances, where it is more convenient to 
refer to the associated density function f, we shall write 
(X,f). In conformity with standard notation 

F(x) = p r(X = x) 

f(x) = 

dx 

Random Numbers : It is considered that there is available 

some formula or mechanism or device v;hich produces numbers 
which are regarded as drawn from a ponulation (X,F) or equi- 
valently (X, f) . Such a number is frequently identified as a 
pseudo-random number in the literature; however, here we shall 
identify it simply as a random number. 

V.’e may say that a set of such random numbers is associ- 
ated w 7 ith the random variable O' , F) or, more briefly, that 
they are associated with the distribution F. 

Senuence of random numbers : As the formula or the generating 

device or mechanism employed produces a set of random numbers, 
these appear serially or sequentially. V-e shall refer to 
these numbers as a sequence of random numbers. 



In referring to random numbers, as defined above, as- 



sociated with specific random variables and their associated 
cumulative distribution functions, we shall need additional 
specific definitions and notation. 



Random T J(0 . 1 ) number : A random number associated with the 

random variable (M,F) where F is U(0,1). 

Normal deviate : A random N(0,1) number which has been obtained 

as a transformation of a random T J(0,1) number. 

Generator: A formula or method for producing a seauence of 

— I, . — 

numbers . 

NORMGEN : The computer program to generate normal deviates 

by the Box -Muller method. 

NOR! IT N : The computer program to generate normal deviates 

by the "sums of uniforms 1 ' method. 

RERQQHN : The computer program employing Marsaglia's method 

to generate exponential random numbers. 

CHI GEN : The computer program to generate chi-square random 



A normal distribution with mean /LL and standard 
deviation U as parameters. 

U(a T b) : A uniform distribution over the interval 

Tn particular, U(0,1) will be used extensively. ■ 

Random N( (J. r fT) numbe r: A random number associated with the 

random variable (N,F) where F is N (// , u) . 




numbers from sums of squares of independent normal deviates, 



with the number of degrees of freedom as an arbitrary- 
parameter . 



i’. 



SECTION 1 




1.1 C THEORY 



For generation of random numbers associated with general 
random variables (X,F) from random U(0,l) numbers, it is nec- 
cessary to establish a relation between the uniform distribu- 
tion and other distributions. Particularly a relation between 
the uniform and normal distributions would be useful, because 
other distributions such as the chi-square, F, and Student’s 
t can be derived from the normal distribution. A general re- 
lation between the uniform distribution and other distribu- 
tions does exist and is called the probability integral trans- 
formation. This relation may be stated as follows: (See [8 

for a proof) 

Theorem : Let X be a continuous random variable 

with distribution function G. Then the random 
variable Y = G(X) has a density function given by 



1 05y?l 



f(y) 



elsewhere ; 



its distribution function is given by 



( 0 yfO 




( y 0 < y «= 1 



1 ySl 



1 



The relation y = G(x) ip called the probability integral 
transformation, and since G is a non-decreasing function it 
follows that (see page 313 of -lfj): 



: r 



G“ 1 (Y) = x] 



Pr 



= G(x) ] = G(x). 



The above provides the basis for generating random num- 
bers associated with an arbitrary distribution by transform- 
ing U(0,1) numbers. 

In the case of the exponential distribution where 
y = G(x). = 1 - e , it is possible to solve explicitly for x, 
that is, x = -log^Cl-y) =0. In this case we say that G has 
been found in closed form , meaning that x has been expressed 
in terms of the logarithm function, one of the functions 
which, like the sine and cosine, we choose to call elementary 
functions . 

For a general distribution function, it is not always 
possible to find the inverse function in closed form, that is, 
it is not always possible to express the inverse in terms of 
so-called elementary functions. 

For the normal distribution, a closed form expression 
for the inverse* function is not available. In such a case 
it may be possible to employ a graphical interpretation. 
Graphical methods, however, have a minimum of decimal place 



accuracy and Lear! to table look-up procedures in high speed 



digital computers. Entering the ordinate of Figure 1 to de- 
termine abscissa values by use of a table, sometimes called 
the "Normal Curve of Crror Table", also leads to table look- 
up methods in the computer. 

d ■ v) 



;/ * 

F igure 1 

The disadvantage of the table look-up method of supplying 
random numbers for a war game is apparent when the memory 
space limitation of a computer is considered. 

Cven though no inverse function may exist in closed 
form, it is possible to construct generating schemes that are 
fast and have acceptable properties. Two such generating 
schemes for the tiormal distribution are given in the next 
two sections. 

Normal Numbers by the Box-Muller Method (NCPMGNN) 

Box and Muller [l] oroposed a method for generating a 
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pa i r of independent random N(0,i) numbers from a nair of ran 
r’on U(0,1) numbers. The transformation (see Appendix A) for 
generating the pairs of N T(0,1) numbers is: 

i, 

X - (-2 log U ) ^ Cos 2 77U 

- 1 “el 2 

t 

X = (-2 log U ) ' Sin 2 7TU 



where 



' T J 

1 ’ 2 



are independent ’2(0,1) and X ,X ace indepen- 



dent N(0,l). 

In view of the four functions (sine, cosine, log^, 
square root) which must be computed in the 3ox-Muller method 
this method cannot be considered fast, although a nair of 
numbers is produced so that computation time must be halved 
for comparison with other methods. The major advantage of 
the Box -Muller method is that it is an analytically exact 
transformation (of pairs of nurbers) from the uniform distri 
button and will faithfully produce a normal distribution 
given a uniform distribution. 



One disadvantage , other than computation time, is that 
the Box-Muller method will reflect whatever bias or contam- 
ination is present in the uniform generator. 

Then written as a subroutine and included in a service 
library, the 3o> -Muller normal generator may be employed in 



two ways : The pair of N(0, 1) numbers may be delivered simul- 

taneously. X in the A -re lister and X in the '-register, or 
1 2 

the subroutine nay rave X until the next call for a random 
number. In the first instance the user would be required to 
save X, for any case where only a single number were desired. 

In general, standard library subroutines are very expen- 
sive in machine time, as they are required to be able to ac- 
cept the general case ...... .they must test for sign and mag- 
nitude in many cases before beginning the computation. In 



special cases where the argument with which the subroutine is 
entered is 0= '=1, or perhaps C = X = 27 T , the tests are not 
germane and precious machine time is wasted. In cases such 
as this the functional subroutines can be rewritten and 
shortened and included as part of the generator subroutine, 
however there is trade-off here, as this greatly increases 
the length of the generator and consequently the library pro- 
gram. The library function routines arc then partially dupli- 
cated within the generator routine when the generator routine 
becomes part of the library. Of the four generators described 
in this thesis, these remarks pertain almost exclusively to 
the V: and ''uller X0T' [GTN. 



1.3 Normal Numbers from luns of Uniforms (NCinSUN) 

\ discussion of the justification for approximation of 
a normal (0,1) random number by the sum of 12 uniform (0,1) 
random numbers appears in Appendix A. The approximation can 
be improved by increasing the number of uniform (0,1) num- 
bers used, but machine time places a restriction on the de- 
sirability of using more numbers. The sums of uniform random 
variables converge more rapidly to a normal distribution than 
the sums of random variables from a non-synmetric distribution. 
The program for generating these numbers is short, in both 
program steps and machine time. There are no library sub- 
routines which must be called in. 

This method is appealing to the intuition as there is 
the factor of "smoothing’' or averaging whatever.* bias or con- 
tamination may be present; in the uniform (0,1) generator. 



1.4 



Exponential 



Narsaglia ' s 



land on Numbers ( EXPO GEN) 

2' method for generating exponential 



ran- 



dom numbers is more complicated in orogramming details than 
either the NOEMGEN or NOTISUU, but lends itself very well to 



the EG 1TEAN language . 

The exnected value of the n of Marsaglia's method (see 
Appendix B) is approximately 1.53, so on the average only 



1.50 U(0,1) numbers are require*'! from which the minimum must 
be selectee , anri 1 . 3 a discriminations (see Appendix B) are 
required to assign a value to n. Assigning a value to n also 
requires an average of 1 . 53 discriminations . Generation of 
an cxnonential number thus requires an average of approxi- 
mately 4.75 U(0,1) numbers. 

This method of producing exponential random numbers is 
sensitive to any bias in the characteristics of the uniform 
generator . 

1.5 Chi-Seuare Random Numbers (ORIGAN) 

Suppose that ten chi-square random numbers with five 
degrees of freedom are desired and that the NCR! TUI! genera- 
tor is us o' 1 to produce the normal numbers. It takes 12 uni- 
form numbers to produce one normal number and five normal 
numbers to produce one chi. -square number wit! « five degrees 
of freedom. Thus we require 50 normal numbers and 600 uni- 
form numbers to produce ten chi-snuare numbers with a rela- 
fciveiy small number of decrees of freedom. Ore has no diffi- 
culty seeing that production of chi-square numbers is expen- 
sive in machine time, however, in war gaming the need for 
chi-s.cuare numbers is likely to be considerably less than 
the need for uniform or normal numbers. 
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The chi-'- ’uc ee * no re tor ’-'rogram is a 5a in relatively 
direct one’ no special difficulties ore encountered. The 
frequency and goodness -of -fit tests do present a special 
problem, as changing the degrees of freedom in effect creates 
a new distribution an' the class intervals and expected fre- 
quencies of occurrence within the class intervals change. 

They do not change in such a fashion that DO-LOCTI may be 
most economically employed. The program becomes suite lengthy 
in machine time if chi-square numbers with several different 
degrees of freedom are tested. 



i .6 Computer Language 



FG IT RAN was the computer language employed for all uro- 
grams written for this thesis. It is realized that FORTRAN 
may not be the best language for all types of war games, al- 
though successful war games have been written in this lan- 



guage. Opeod and memory space limitations may dictate the 
use of assembly routines as opposed to FORTRAN - type comnilers. 
However, within a FA '"’RAN program, strings of symbolic machine 
language may be inserted where the legal FORTRAN statements 
are not appropriate. In general, the particular character 
of the simulation will suggest ;he most convenient language. 

For use in a war game, the four random number generators 
in this thesis must be translated into the language in which 



tho game is written. One wQr r ’ of caution is in order at this 
point. The nethod used to generate the U(0,l) numbers de- 
pends on the computer having a 43 -bit word length, and the 
method must be revised for use on a machine with a different 
word length. 

The Control hate Corporation 1604 stored program, gen- 
eral purpose digital computer with a storage capacity of 
32,763 4° -bit words was available at the Naval Postgraduate 
School for this project. 
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2.1 oner" 1 Fa c’ :g ro* mH 

The ourrose of this rection is to discuss the following 
question: ' her, a generator has produced a sequence of ran- 

dom numbers and - t is desired to associate then with a parti- 
cular distribution function F, what should be the criteria 
for deciding on toe merits of toe generator? 

For illustrative purposes the X(0, 1) distribution will 
freouentlv be referred to in this section, although most of 
the remarks will pertain to the generation of any sequence 
of random numbers. 

The criteria will to some extent depend on the employ- 
ment of the generator. For example, statistical sampling ex- 
periments for estimating the parameters of mixed populations 
require vcrv strict accenta.nce criteria for the random number 
generator. The employment considered in this paper will be 
nrimaridy in the war game area. 

Tf the F with which the sequence of random numbers is 
to he associated is >’(0 f l), it is obvious that a mean near 
xero and a variance rear one must be demanded, but is this 
sufficient' 1 The answer is enohe.t 1 ca 1 ly no, as Figure 2 

j j 



i lluctrates . 




Figure 2 

e Tensity function graphed in Figure 2 could well have r 
near, of ~ero and a variance of one, but the skew to the left 
elinina**es it from considerat ion . 

A teat of the sequence of n vibers by computation of the 



first three central nonents gives the criteria for the mean, 



variance and skew . Is this suffiecient? The answer is defi- 
nitely no again. Two density Sanctions nay have the sane 
values for the first three central moments and yet have quite 
different values for the fourt’ and higher central nonents. 



urtosis 1 s 






: qeakedness' T or '"flatness 1 



of a 



en- 



sitv function. \ ro .igh measure o e 
a a’ka'- 1 fourt 1 '. central noser t . A 
tively high res’ such as the curve 
lonto’-urtie, T . i 1 c *-he curve in Fi 
"o-eped is called nletykurfic. \ d 
e ,r cessinf 1^ nan 1 '®' 1 nor excessively 
tic' , fo” example the cornel d i s v r 



v uetcs is is the stnnd- 
d- stribution having a re la- 
in Figure 3a is called 
g ire 3b which is flat- 
istribution which is neither 
flat is called 'hpsolur- 
ibution see 111 . 
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Flatykurtic 



Leptokurtic 
Figure 2a Figure 3b 

o o 

Now we have established a requirement for at least com- 
puting the first four central moments to provide data for de- 
ciding if a sequence of numbers can be regarded as a sample 
from a desired distribution, F. 

At this point the experimenter would do well to apply a 
frequency test to the sequence of numbers by dividing the 
theoretical range of the numbers into class intervals and 
counting the numbers which fall into each class interval. 
Since the theoretical frequency of occurrence of numbers in 
each class interval is known, a goodness-of-f it test like the 
chi-square test may be applied to test for a significant dif- 
ference between sample frequencies and observed frequencies 
in the class intervals. 

In the case of the N(0,l) distribution, any sample con- 
sidered to have been drawn at random should have an approxi- 
mately equal number of positive numbers and negative numbers. 
A count of positive signs thus provides the data for another 



r: 



tes 


t. Tif occur 


re nee o f par i t ive 


signs for any 


distribution 


has 


the binomial 


distribution wit’i 


p - h, if the 


median of 


the 


distribution 


is at aero. dee 


feet ion 2.4. 





The tests dis cur> see 1 up to this point were concerned with 
the distribution associated with the generated sequence of 
random numbers. The order in which the numbers in the se- 
quence are generated is also of orimary importance. Tor ex- 
ample, in the general case, a sample of say 10,000 numbers 
could be generated which would pass all six of the tests for 
a given distribution and yet have been produced as a nonotonic 
sequence. Olearly the war gamer -would not ordinarily be in- 
terested in tliis sample. 



how then do we test for randomness, in contrast to test- 
ing for the characteristics of the distribution? Since we do 
rot want the numbers in the sequence to have a predictable 
relationship t:o numbers following later in the sequence, a 
test for the serial correlation of the numbers (with various 
la g s ) is d c o irablc . 

The numbers may show a low serial correlation and yet 
co-' a ‘ excess ivo’ " - runs of either oositive or nesdr'v 1 

»' '• *s. 

lumbers. Tor the war gamer this could have disastorous ef- 
fects, as. his rx-x cted values, could be based on a series of 
’?c : ion -.oints • -V<ro Ms random number generator is rroduc- 



1: 



in. - ; a run of 
The over-all 
be player 1 a 



all positive numbers , or all negative numbers, 
effort of this situation is that the game must 
prohibitively large number of tines to obtain 



results that have a meaningful average. 



Table 1. 
a ary tests 
proper tl es 



bovo general testing information is summarized in 
Tt is to be emphasise.- 3 that these tests are neces- 
, but are not sufficient to guarantee the desired 
, even when administered collectively. 
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n explanation of 


"he meaning given to , the sample site . 
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’ace.; the odd numbers between 1 and 2 


These numbers are 
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transformed -r.t-j i i i ^rval (0,1). The 7 J(0,1) numbers r re 
transformed in one -to -one correspondence to numbers of another 
hist ribution, so a cycle in the generation of the 0(0,1) num- 
bers implies a cycle of the sane period for the transformer’ 
numbers . sample rice of n - IOC means the c irs t 100 icun- 
bers of the cycle arm’ a sample sice of n = 500 means the first 
500 numbers of the cycle. In the second case the first 100 
numbers a re refe^ed as the first part of the sample of 500 
numbers. In of he** words, all samples start with the first 
number of the cycle produced by the uniform generator. 

Oince some °, 0 Q 0 billion numbers are produced before re- 
peating, an euperimenfer could enter this cycle at any ,3 c- 
sire'- 1 point and start qeneratinq numbers which suit his pur- 
poses. The tests in this thesis extend only to the firs*: 
10,000 numbers of the cycle because of the restrictions- of 
available machine time. The generation of 10,000 normal num- 
bers is ~<uite rani.d but the batterv of tests aoelied dras- 
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tically increases computer time required. 
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The moment tes-r. are constructed from the definitions 
of **he first four central sample moments. /-II summations 
are taber on 1 , fo 1 * •- 1,2, ,r. 
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n, through n, are conouted from venerated sample. values and 
are consistent estimates of the theoretical moments. 

Ml moments exist for the normal, exponential and chi- 
square distributions, and the theoretical values may be ob- 
tained from the moment generating functions. For the N(C,l) 
distribution the odd moments vanish and the second and fourth 
moments have a value of one and three, respectively. Compar- 
isons of the first four sample moments with their theoretical 
values re made for all four random number generators in this 
nap or . 
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where m , m and, n are the central moments defined above, 



lq and Q> are called the coefficients of skew and hurtosis, 
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tV di 5fe'*cnce 1: s *n the choice of too class intervals. T r 
the first case the -bscissa is divided into 14 class interval 
as shown in Figure 4, Twelve of the intervals are of equal 
leapt’ an.’ the numbers creator than three arc combine'' 1 into 
a simla class and the numbers less than minus three are con- 
bin'*'- 1 into a single c 1 acs . 
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Tor the second test *-hc class intervals are assigned 
equal probability of occurrence of sample values, rathe" than 
O'-ual len.qf'h -vibd-iv- siens . The class intervals are chosen as 
indicated in ’"inure 5. 
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"or the exponential and chi-square distribution goodness 
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erial Correlation Test 

'"he serial correlation coefficient ir a measure of the 

relationship o c m .iba„rs in a sequence to numbers occurring 

later in “he sequence. For example, given r sequence of n 

numbers , ’> ) the relationship of the first nurr.be 3 

12 n 

So H'i r r, > the second to the fourth, the thir^ to the 

fi ft*', r " r ^ so on ir '-ailed serial correlation, lag tv:o, and 
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-i ’ nayeti'*e with no runs £ tea ter then one in length above . r 
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•’or. , <_r. t’ e r ~ ‘ ■ '■ul * re ‘-oo yoo ’ . The greatest tv tuber of 



.nr r 'O 
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;*>■' . Tn f 'h’eir setinle the runs nay be counter’ enr? 
connote J ” ? f-u ;-n r t-heor *' ■' o' - 1 - >p 3 u ftr for Jr to minin' - * siop.if i- 
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oT Slori b jlj, t’ r number of positive an' 1 n ega 
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'"ince it ir ^es ired to keep t;be expected number of nr 
in any class interval greater than some number, say 10, for 
- goo-Y;Cf s-o*-f i "ert, the expected numbers of runs greater 
than ^ certain length me 3’ be grouoeO info a single class. 

The (5-^^ •’ numbr .• j r positive anf negative runs of length 
renter than or equal to some number k is r.^roxima tely : 
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Let he refine'’ bv the following equation: 
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r at least one run o c length? j = (X 
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r suits o c the t~‘sts performed on sequences of n.tm- 
!-ers venerated by '’.?T~7N, y-pyOUh, ’rjOGCN and CHIG5K are 
given ‘ ’’’abler 2-*'. ~t rbould be noted that ^os Corr. and 
No 3 "orr. in ”>! !es 2 - r refer to the largest positive and 
largest negative values of serial correlation obtained in 
test, r o r lap me thr nqh 20 . 

7ur the normal venerators in Tables 2 and 3, C’IIGOO .1 
refers to a fre * 1 lencv test with e ual class intervals and. 

de'?"es of f" ci, an’ 0710 *' n 2 refers to a frequency test 
with equal exp *> 'ted frequencies in class intervals and nine 
deprees of ^reedoi . The theoretical values of chi-square are 
the valuer above which t- r test results are significant at 
the ,05 level. 7CTt77N and M 0 'T-°Uu both oar see’ the goodness - 
of-~it *" ■ r or run-' *r ra-nl^r of size 1C00 and greater at 
the ,05 level c r significance as indicated in Table °, arid 

.os 

as define ’ 'v Occ.t ion 2 . 

-s h 00 7’ Cf.i is to pass the frequency tests for sample 
size o c 5000 and lo^qn at either the .05 or .01 significanc 
levels, the frequency count for the ten class intervals for 
cannier of rim 1000, 5000 and 10,000 are compare ' 1 in Table 6 



there were no runs longer* than K , where 1’,,^ ~ 3. 3 Clog n 1 



intervals are 



.‘.lthoug 1 * the perceritag ( er of numb.- as in class 
similar, the fit becomes poorer as the sample size increases 
due to th.e square difference in the numerator of the chi- 
s cur re statistic. For example, 120 x 120 ~ 14,400 is 100 
timer as qreat as 144, while 10,000 (the denominator) is only 
ten times creator than 1000. Meo it is noted that a sample 
of size 100 fares very badly in most tests performed on se- 
quences of numbers produced by T'.rOGFN. The magnitude of ex- 
ponential numbers is very sensitive to the magnitude of the 
uniform numbers from which t’^ey are transformed. This is true 
of harsaglia’s method £oj and a more we 11 -known method, i.e., 
taking the logarithm of U ( 0 , 1 ) numbers. 

The last: col rra of Table 5 (C’lIOFi!) gives the values of 
th.e chi-s.qjare goo<’nr ss-of-f it statistic using 10 classes. 

For nine degrees of freedom th.e . n 5th and . on th percentile 



points of the c^i- square distribution are 16.0 and 21.7 re- 



spectively. Thus, significant results at the 57 level -were 
obtained for r * °^0 and U - ?. ,4,5 and 15. In order to in- 
vestigate this further a chi-square analysis was done for 
sarnies of size 50'd and 1000 (Table 7) . It should be noted 
that fit does not necessarily improve with sample size. For 
all degrees of freedom and sample sizes considered, CHOFF 
passes the test of fit at the .01 level of significance. 
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fortV * e-%, lc»*e the properties. of the O’lIGa!’ numbers , 
the fm rsf feu * ior (. ids were computed but not tabled for para- 
plot cf pace l'"' 0 , 500, 10^0 and 5000 with five degrees of 
freedom. The t’eoretfcal values of the first four moments 
ere ,1'*, 4-0 and dO, respectively, as commuted from the 
moment generating function. The values of the first and 
second •lomentr were reasonably close to the theoretical val- 
ues, but the values of the third moments were quite low for 
all sample si.es. The value of the fourth moment was low for 
a sample size o f 50'"', indicating a possible ’’flat' curve for 
the >ns ,? tv. 

"one of the raw data test results have been included in 
"'dies 2-~, as they are likely to be of interest to experi- 
menters. Table ° is a summary of results of tests for n = 1000 



:or 



r>* O’-tvi 



an'’ n - d, 000 for all other generators. The terms 



G ’'O', T.n i nd in Table ° are used in an intuitive and 

relative sense rather than as any absolute standard. GCd 
generally implies that the theoretical and test values are not 
sign if leant lv different at the .05 level and no undesirable 
characteristic? are apparent. TOGO, implies significant dif- 
ferences at the ,01 level or sufficient deviations from theor- 
etic'! v lacs to take *:he generator suspect for some applica- 
tions. FV-.IT implies a condition somewhere in between the two 
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The teste have Semens t rated that random numbers nay be 
generated that behave when tested like numbers from a normal 
exponential or chi- square distribution. It is apparent, how 
ever, that the behavior of these numbers depends heavily on 
the multiplicative scheme used to generate the T J(C,1) random 
numbers from which they were transformed . home of the gen- 
erating schemes discussed in this thesis are theoretically 
exact so that t u e tests described in effect provide addition 
al tests of the performance of the random. U(0,l) generator 
that is an integral ; art of NCTfGdN, NOSHSUM, T'PCGTN and 
C T, IGd?7. for generation and testing of random 

number^ . 

In general, the test results of Section ? indicate that 
V" T'.G'Id r'V ITTTd .n jst be used with caution for verv small 



samoies and N 



^ v must be used with caution for large sam- 



ples. The chi-square generator was obtained by two transfer 
ma •"buns -pd should be a re fully examined for the sample sice 
and denr-ecs c r freedom required, as the results are not con- 
sistent throw .houfc. further testing is required for the 
latter and will require rather extensive machine time if 
large samnle <• and <?pveral. numbers of degrees of freedom 



°1 



re cor..' i 



• ’--.i 

nm numbers r 
these - \tnberr 
^ or'' rotor has 
in 7; a aero in 



r . , -n; ci~>cnt) which uses a ;tent many rar- 

oul •’ ’ ic^e the capability o p testing samples of 
at des'rod intervals to ensure that the number 
not tv 1cm into some sort of trap such as find- 
a Tultiolicativo scheme. This, of course, is 



merely an extension of the general rule of not trusting the 
machine output answers before submitting the morel to a test 
prog "an with hnovre answers. In any extensive program the 
arithmetic rc v .enc“ may ' e giving answers to questions that 
were really a she*’, n ' not to questions the programmer thought 
h-' as' od. 

.\r an incidental observation, it is believed that the 
logical sequence of learning program techniques is to start 
with the machine language and the fundamental machine instruc- 
tions, ovre^' 1 to -a assembler anr’ finally to a compiler. 

..van t ! 'ou; h a non"’" 5 lor language such as ^CPT’V.M is advertised 
a- "no computer <? r-erience necessary", the programmer who 
floes noi unhe ."'tan ’ » u c basic machine instructions will fin / 1 



•V-1 urging M r nron'c- a nearly inpossible tas’' . This is 
apparently not a ♦'ter of e-neral agreement. 

lince random numbers have- been generated by several 
me. ' r.tVv *-ha*o a ‘ bmetic schemes, including generation by 



.1 J\ ..5 1 «.i: ' ■ 1 .t led devices i , i** would seen reason- 

able to make a c- n. .ricon, b-zo or three of the "best’' pub- 
lished table? of random numbers could be typed onto cards and 
inserte J into machine memory for testing by the same battery 
of tests that r ri h, MCBMBUM, BbIG'dN and uTOfrih were sub- 
jected to. It is real 'red that typing 10,000 IBM cards is 



me. .. ! * as 



T_-p s on o c.ci c o °* 



ren- 



^ '’ecks of car 'c, mar be available 



The iterative, multiplicative algorithm used to gener- 
ate the ’ r ( A , 0 ce nenc" delivers a long cycle of numbers vritV 
out reporting i nd* vidua. 1 numbers. Blochs of these number’s 
within th"' cycle re suite lihely to have undesir cable char- 
acteristics. The cycl"' can be entered at any point by in- 
sert* nr "ho Vs hr'' ’'striding number 1 ’ in the algorithm. If 
small bleed's of numbers were tested, those blocks showing 
no ~ r f’iir: ctari- ! : les could be dropped from the cycle by test- 
ing for the critical priming numbers. There is a definite 
d,- rimer involved this procedure . The drooping of short 

ilncl? of numbers '-u.ng poor characteristics may improve 
c v . 1 1 samples and ’a the characteristics of large sam- 

ples. 7’ is matter has apparently not been discussed in the 
literature . 
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ables by performing discriminations on the relative mag- 
nitude? of uniform (0,1) random variables. 
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Anals of Mathematical Statistics 
Volume 32, Number 3, September, 1°61 



Suggestion for expressing a distribution function as a 
rn Sure of suitably chosen distribution functions .... 
choosing a distribution function which is close to the 
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correction only infrequently. 
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The purpose of this* appendix is to discuss in more de- 
tail the two methods (NOT 'SIN an d NBTtSUN) of generating 
random *’(0,1) numbers introduced in Section 2. 



it Xf 



•puner 1 , as has been rone throughout this thesis, 



that r. method to generate uniform random numbers is available, 
0 transformation of uniform random numbers to normal random 
numbers will now be considered . Of the several existing 



methods to accomplish this, the method suggested by Box and 

'duller 1 1, was selected for accuracy and integrity in the 

tails of the distribution, as war gamers and others can be 

-ui to interested in "rare events . 

NCnSSN: Let 11 .U be independent random variables from 

1 ’ 2 

the uniform density function on the interval (0,1). Consider 
the random variables: 



(-2 log 7 J ) Cos 2 T r U 
e 1 2 



( 1 ) 



X 



(-2 log U y Sin 2T7’J 
"el 2 



Then 



(X, ,X ) are a pair of independent random variables from the 
i 5 2 

same normal distribution. Notice that the radical is always 



positive since 1. 



Justification: 



Solving (1) for the inverse relation- 



s i. 



w* .tV>ut dif ficulty : 



- e.o I- xh /2 

I - X /. 



- ( 1 / 2 7 T) arctan (X /X ) 

2 1 



tT* 1 T i • 

to iaco:: un oft-. 



, , )w ith respect to(X, ,X^j is found Co be 



(1 /2TT)w> [-(XV XV)/2 



1 



and the joint density function of 



bh, T 0 is 1, so that 1 he jo* nt density function of X, ,X is: 

a/2 7 r )eno [-<X* + X^)/?.' 

(1/2H V<o LxV? . ( 1 / 27 T) 2 expf“'C /2 



c ( v ‘ 

■ V v 2 



f(x )• f<X ) 



•'rtrbl ' sMn f lat X, and X are normal ( 0,0 and independent. 

>:oress ’ on ( 1 ) is analytically exact and the character- 
ise - ] cs of the r.urbrs produced are distorted from the ideal 
only ''V t 1 e pseudo -andom uniform numbers used as input and 
by th ~ ia chine :a 1 a- ' • a t ' ons , inc lu<* ins the square root , lop , 
sine and cor.irv •f'uti-op . T b:nected accuracy is 13 decimal 
djq-* r <5 ’ n the CT >0 If 04 70 ' TT?,.XN routines. It is hoped and ex- 
pected chat connu*’''tional round-off errors will be random. 

'V l wc l!T : V ' ie users may not be interested in rare 

events but ir Timbers around the mean we next select for 






to? fir met’ o' 1 a red on the Central Limit Theorem. The 

■distribution of on e of independent random variables from 
symmetric distributions approaches the normal distribution 
guite rapidly. It is noted that 50 numbers, each composed 
of the sum of five random dibits, gives a strikingly normal 
appearance when plotted in a histogram. 

Cince the variance of the U(0,1) distribution is 1/12 
and variances of uniform random variables are additive under 
convolution it is convenient to select 12 as the number of 
uniform random variables whose sum will approximate a normal 
random variable, beans of uniform variables are also addi- 
tive so that it remains to subtract the constant six from the 
sums o f 12 independent 11(0,1) random variables to approximate 
a normal (0,1) distribution. 

This process truncates the distribution so that -6 = x = 6. 
1 low much is lost by this truncation can be estimated by the 
fact that normal tables with 4-place significance have posi- 
tive probability within the limits -4.2.6= x =4.26. 

There is a cast 1 f or the truncated normal distribution 
in that a model is useful when it can be realistically ap- 
plied to a real world situation. Consider the height of men 
of a civen race as a random variable, normally distributed. 

"or use the distribution is truncated in that the smallest 

' r 



' : ht tr c t. ’ 1 ••.re; fc-r than one foot one’ t u e greatest 

height lest than ten feet. Definite bonne’s can be set on 
many errors of nea ummont which are considered normally ’is- 
tributed. hnpl : cations of the theoretical normal distribu- 
ted are so frequently truncated in the real world that it 
seems inconsistent to reject a manually generated sample from 
a normal population if the sample exhibits the properties of 
the theoretical distribution between the limits of usable 
values . 



APPENDIX B 



EXPONEN TIAL RANDO M NUMBE RS 

The exponential distribution figures prominently in 
particle or radiation studies, reliability, life testing and 
like areas which can be of interest to persons using high 
speed computers and monte carlo techniques. A power series 
expansion to compute the logarithm of a uniform random vari- 
able may be used, to generate exponential random variables, 
but a more rapid method is desired. 

0. Marsaglia of the Boeing Scientific Research Labora- 
tories offers a simoler .device for oroducing exponential 
random variables by performing discriminations on the rela- 
tive magnitudes of uniform (0,1) random variables. Bee '2 . 

The idea is to choose the minimum of a random number of 
uniform random variables, then add a random integer say. 



let n and m 


be random 


integers 


according 


to the following 


schedule : 
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■ ■ :-j ,M , 

i ^ 



,’J ) 

n 



h<- s the oxoopo i al 3 \ <?tribution. 
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ien 



b i J ; r i r. ' : I t n be a re. , fori variable tahirv; 

•" 1,V, r'f' ~ robrbili.ties p ,p , 

I 2 J 

- - oi (h , j , >1 : 

1 2 n 

tb -1 flrt ibati o c y is, for 0 = 9 = 1, 



l-!y<l-«; - -y !-*}**- 



;v -tic-lb •, if C 1/Cc-l) - .5*12767 



. . . one 
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c/ - • > 



ic r . 



*“* «n ' ■ 



valuer ^,1,2, ... 



re ( 1 -o ) , 0 ^ = 1 1 

T-> e o r •■".! : Tf c - l/'c-l) nnfi if “he lrn^on var- 

iable a tehes *.v 1 ier 1,2,2, ... with probc.bili- 
5 ■ r c , c • , • , ... <rr ,! if, iiv’epnnf ently, 

- ran* 'em wr - ' ’ 

;oK i'ri •? l/'co), l/(ce”), l/(ce‘ ' , ..., 
t'vn t v e rarvh n y* 'able 
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’ <- th<? ev o^-er** 4 ^ * 1 i r t ribut ion 
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~roef : a ~r • ter o* verification if v = 1 
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ait enpix c 



:hi square pantxt: _ muiib prs 

Let ... , x.^. be y independent random variables 

normally distributed, each with zero mean and unit variance. 
. 2 

Let 2d equal the sum of the squares of these \J random 

2 2 
variables, i.e., Y‘ = x* + 'X t- . . . + . X is also 

a random variable because it is a function of random vari- 
ables. The density function of the chi-square random vari- 
able is written as 



1X1 



*i = 



? inn 



. ife~) .x 

(XT X 



z 

X ^ 0 






otnerwise 



All x's must have zero mean and unit variance. If they are 
independent normal ( // , .X) numbers, the quantities 



A- 



- * i.2, ,1/ 



do have a normal distribution, each with zero mean and unit 
variance. Hence, --urn; iny on i, for i - 1,2, , 'j 



X' 



= eh : --i\y-XV<d 



has a chi-square distribution with 'j decrees of freedom. 

An example of such a random variable is as follows: 

A missile is sent toward a target. If the range and de- 
flection errors of the missile are independent normally 
distributed random variables each with zero mean and equal 
variances, the squared distance between the explosion and 
the target, divided by the variance, has a chi-square dis- 
tribution with two degrees of freedom. This can be easily 
extended to three dimensions. 

A subroutine available to a war gamer to supply inde- 

2 

nendent numbers with arbitrary degrees of freedom thus 



can Ti ave immediate aonlication. 



A3 1 :;r i r ' 

^ > a T 'C* 

__ -V W. X . > 

Following is the Por-Nuller method of generating pair? 
of N(n,l) numbers from ’-'airs of 0(0,1) numbers, programmed 
in r^e Control Bata Corporation 1^04 FORTRAN language. 

The program r r ’ written merely computes 100 numbers 
and stores them in an array called ''TAB”. K is the sample 
size and V must a lx. -ays. be one-half the value of N. The di- 
mension of TAP must be at least the value of N. The numbers 
that the program stores in TAB are in floating point format. 



"T 'TNOIQN I 7 J ( 2) , TAB ( 100)" 

!tASK=4000C0rtOOO0Onn003 

'.'12 -?..*? . 

UNIT--- 77 .**22 . 

I UN TO JNTT 
T’(2 )-l 

M=10O 

T ?0 N= 1 , 50 

FVI2(1) , INI 1(2) ,L r '’Al(IU) ,NUI(IUNIT) ,SCL(:iA r K) ,STA2(IU) , 
NIJT(IUNIT) , 
c.clOta::) . 

0TA1 (T T J) . 

T J1~TU( 1) 

’.Tl-Ul /2 ,' v v 47 
U2-TU(2) 
f 2 -3 7 2/2 7 

'uao=i ./vi) 

a -mo 7- toA-tto 

' T vt O OTM^/ A -n^ \ 

JU* /_ - )ih \ . t. x i , 

n ' ry _ r* '> n r? / ^ nmN 

' > • “ » V- -> \ V J 

v T ''/“I'* 

;i= roT-ocop 

t — V' v. X > JL *. . t 

T'-OtV' ' 



f *7 
• / 



?o 



L -1 

tab(D= :2 

TAB(L)= :i 
CCCITIUUK 



The TYl lf n 4 nrogram for generating N(C,l) numbers from 
the sums of 12 uniform numbers also stores M numbers in ''TAB". 



r R' s RR A ' f NOR 
or KINS ION TAB (100) 

N=100 

MA S K= 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 B 
UNIT=2.**23+3. 

IUNIT=UNIT 
VJ - 1 

or 702 t=i,:: 

BTC 0=0 

00 701 ”=1,12 

LDA( TU) 'UI ( IUNIT) , SCLC 'ASIO , STA( IU) . 

:iui=iu 

aui=Kiui/2.— 47 

701 3T0?«.= OTOK+ ' : IU 1 

rr» ▼') e> rv> /- r*\ 

Jiv A=.'iv-V- • * 

702 T/.B(I) = ST0K 

The chi-square generator has tv?o parameters, the size 
of sample N and degrees of freedom NU. 

m ^ f rsi r -r r\ T 

- v*. ^ Vt i. . - 1 0“ 4.S 

^riSNSlON TAB (IOC) 

MASK=40000000000000003 
'X\ 7 T1=2 .—23 ‘r ? . 

NU= 5 
N= 100 
IUNIT=UNIT 
I1J=1 

o: °C3 i=i,k 
ITOTlsO 

OQ 302 TT=1 UU 
"TO 3-0 

oc 301 !'=1,12 

LPA ( TU) Til ( TUN IT) , 3CL(: CASK) , C TA( IU) . 

kiui=iu 



4.2 



rH 



i ; ; i - - 1 



T' 

n 1 _ 



J. ~ 1 



- 0? 



x th ,T \? 7 



TX/lP/TX! 



/V - v j i. s 



rt '.e exponential ;encrator ’las th® ainple parameter N. 



JL ^ J . v- 



n i; ' xt j i c n t.\" ( ’. o°) ,tx! C? o) 
:-'A'?TC=20000000000000003 
! XT<2=XX0OX00000000F 

▼ o y j r*o 

jNIT= 2 . ** 23 f 0 . 

t:=ioo 

•^T t\T Ym t y* y 

L * » J 1 — J s 

c r -i . / ( ' : - 1 p 
ir=i 

X 1/3 L-l.N 



LDA(TU) , 

:rJ T (iuNi't) 

r>rr-. a / TT v N 

*„i ! xU ; 

.X' 11? 

'xcxxio 

» r-^ /» a a r ^ 

N * ‘ ■ *• ' 



• ””=1 

i?/- xt -!*:”■ 

T 



T 



N 3 r- *» o r pc 

' •' ■ > 1 5 



•"* ~ “ nr*/ 

(r.i) 

■'"lO-.'NTT' 

2 3 LC "/:’ 2 N 
"TprO 
' V (10 

* -vp '' - » ?~\ 



7.--C 
°T A ( 0 

3= 1 . /cx 

H=0- 



137 


3=3/7 

T T __ T T , 




Tn/u T T v 1 o 

ir Ji-A, 1 ' 


136 


’ *"» f T f - 1 
•J l— ... ‘ i # 




CC "C l ' 5 7 


13° 


L?A(: i:: ) 




''lATT / *5 or\7 



:ni?(i) 

13° 1 L3A(TU> 

AUDium) 

ccr ' r AH'?) 

?TA( IU) 

A 33 ( 11 ) 
a^C 7 .cr) 

FAA(MA?:0 

?TA 2 or;?: r ' 

1 ^ fj TOT’-O / ^ r> \ 

.i. W J . -L i \.z_ \ * / 

?LJ( 1°°1) 

oc 141 r-2, !':•■:• 

T 7 (T’Tr(i)-T , ::r(: 0 ) 141,141 , 142 
142 T'r'{l)=??' r A') 

141 ccrTTMj': 

14 ? ” 7 . 3 ( l) =? 

v ' o C3C K O' m.mfor test orogrnnis written for use in 
this thesis are available in 7 ' it A7 language in the Computer 
Center of the ”, 7avnl "ostgra^uate School, Monterey, Cali- 
fornia . 
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